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ABSTRACT 



Context. Low-mass stars are element factories that efficiently release their products in the final stages of their evolution by means 
of stellar winds. Since they are large in number, they contribute significantly to the cosmic matter cycle. To assess this contribution 
quantitatively, it is crucial to obtain a detailed picture of the stellar interior, particularly with regard to nucleosynthesis and mixing 
mechanisms. 

Aims. We seek to benchmark stellar evolutionary models of low-mass stars. In particular, we measure the surface abundance of '^C 
in thermally pulsing AGB stars with well-known mass and metallicity, which can be used to infer information about the onset and 
efficiency of the third dredge-up. 

Methods. We recorded high-resolution near-infrared spectra of AGB stars in the LMC cluster NGC 1978. The sample comprised both 
oxygen-rich and carbon-rich stars, and is well-constrained in terms of the stellar mass, metallicity, and age. We derived the C/O and 
'^C/'^C ratio from the target spectra by a comparison to synthetic spectra. Then, we compared the outcomes of stellar evolutionary 
models with our measurements. 

Results. The M stars in NGC 1978 show values of C/O and '^C/'^C that can best be explained with moderate extra-mixing on the RGB 
coupled to a moderate oxygen enhancement in the chemical composition. These oxygen-rich stars do not seem to have undergone 
third dredge-up episodes (yet). The C stars show carbon-to-oxygen and carbon isotopic ratios consistent with the occurrence of the 
third dredge-up. We did not find S stars in this cluster. None of the theoretical schemes that we considered was able to reproduce the 
observations appropriately. Instead, we discuss some non-standard scenarios to explain the puzzling abundance pattern in NGC 1978. 

Key words, stars: abundances - stars: AGB and post-AGB - stars: evolution 



1. Introduction 

Modelling the final phases in the evolution of a low-mass star is 
a demanding task. The stellar interior is a site of rich nucleosyn- 
. !^ thesis, particularly when the star evolves on the asymptotic gi- 
, ant branch (AGB). The freshly synthesised elements are carried 
^ to the outer layers of th e star by me ans of the third dredge-up 
(TDU, for a review see iBussoet all [1999) . The onset and effi- 
ciency of this mixing mechanism depends on the mass and the 
metallicity of the star. Various stellar evolution models devel- 
oped in the recent past (e. g. j Straniero et al. 1997; Herwig 2000; 



Stancliffe et al. 2004; S traniero et al.l2006tlKarakas & Lattanzid 
2007 ) agree on the theoretical picture of the TDU, but details are 



still subject to discussions. It is thus important to test the mod- 
els against observations and to derive constraints to improve our 
understanding of the involved phenomena and related problems. 
AGB stars are tracers of intermediate-age stellar populations, 
which can only be modelled accurately when the evolution of 
the constituents is known. Moreover, evolved low-mass stars un- 
dergo strong mass loss in the late stages of their evolution. They 
enrich the interstellar medium with the products of nuclear burn- 
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ing and thus, as they are large in number, play a significant role 
in the cosmic matter cycle. 

The primary indicator of TDU is a carbon surface enhance- 
ment. An originally oxygen-rich star can be transformed into a 
carbon star by an efficient dredge-up of '^C to the surface. The 
criterion to distinguish between an oxygen-rich and a carbon- 
rich chemistry is the number ratio of carbon to oxygen atoms, 
C/O. During the thermally pulsing (TP) AGB phase, while the 
C/O ratio is constantly rising as a consequence of the TDU, 
there are mechanisms that may counteract the increase in the 
carbon isotopic ratio '-C/'-'C. The radiative gap between the 
convective envelope and the hydrogen burn ing shell can be 
bridged by a slow mix ing mechanism (cf. 'Wass erburg et al.l 
1995; No lle tt et all2003l) . In this way, '"C and isotopes of other 
elements are fed back to nuclear processing. The physical origin 
of this phenomenon is still not known, which is another argu- 
ment to establish observational data in order to benchmark theo- 
retical predictions. 

Measurements of the carbon abundance and the carbon 
isoto pic ratio have been done for a few bright field AGB 
stars jLambertet al.ll 19861; [Harris et al.ll 19871; ISniith & Lambert! 
Il990h . The direct comparison of the values inferred from field 
star observations to evolutionary models is complicated by the 
fact that luminosity and mass of those targets can be determined 
only inaccurately, while both quantities are crucial input pa- 
rameters for the models. A strategy to circumvent this prob- 
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lem is to observe AGB stars in globular clusters (GC). They 
provide a rather homogeneous sample with respect to distance, 
age, mass, and metallicity. Admittedly, this simplistic picture is 
slowly disintegrating. There is evidence that globular clusters 
often harbour more than one population with abundance vari- 
ations among the constituents ([Mackev & Brobv Nielsenll2007t 
iPiotto et al.ll2007l : lRenzinil l 2008l) . The parameters like mass and 
metallicity are, however, still much better constrained than for 
field stars. For our purpose, i. e. investigating the abundance vari- 
ations due to the TDU, the old Milky Way globular clusters are 
not well suited. The AGB stars of the current generation have an 
envelope mass too low for TDU to occur, while more massive 
stars have evolved beyond the AGB phase. The intermediate-age 
globular clusters in the Large Magellanic Cloud (LMC) serve 
our needs better. We started with an investigation of the AGB 
stars in NGC 1846 (Lebzelter et al. 2008, henceforth paper I) and 
indeed demonstrated the abundance variation along the AGB due 
to the TDU. 

In this paper, we pursue this idea and study the AGB stars in 
the LMC globular cluster NGC 1978, deriving C/O and '^c/i^c 
ratios. We give an overview of previous studies of this cluster 
and our observations in Sect. |2l Details about the data analy- 
sis are given in Sect. [3] In Sect. |3.4| we describe the evolutionary 
models we use for a comparison to our observational results. The 
results follow in Sect. H) We discuss our findings in Sect. |5] giv- 
ing scenarios on how to interpret the puzzling case of NGC 1978 
before we conclude in Sect. |6] 



tidal interaction with the host gala xy seems to be a mo re likely 
explanation for the high eUipticity dvan den Berghl2008l and ref- 
erences therein). 

iMucciarelli et al.l (l2007ah derived the cluster age, distance 
modulus, and turn-off mass by fitting isochrones to the clus- 
ter colour-magnitude diagram (CMD). The resulting values are 
T = 1.9+0.1 Gyr, (m-M)o = 18.50 and Mto = 1 .49 Mq, respec- 
tively, for the best fit found when using isochrones from the Pisa 
Evolutionary Library (PEL. ICastellani et al.l2003l ). When adopt- 
ing other isochrones with different overshooting prescriptions, 
the parameters vary so that one could also derive a somewhat 
lower turn-off mass (Mto = 1 -44 Mq), a lower distance modulus 
((m - M)o = 18.38), and a cluster age from 1 .7 to 3.2 Gyr. 

The various radial velocity measurements quoted in the lit- 
erature widel y agree within the error bars, e. g. lOlszewski et al.l 
61991,) and .Schommer et al.' ('1992') derived (vj = H-292 + 
1.4kms"' whereas Ferraro et al. (2006) determined a mean he- 
liocentric velocity of NGC 1978 of (vv) = +293.1 ± 0.9 km s"' 
and with a velocity dispersion cr = 3.1 km s"'. 

NGC 1978 har bours a number of red giant stars 
("Lloyd Evans' '1980) some of which are known to be car- 
bon stars (Frogel et al. 1990) . The easiest explanation for such 
an occurrence is that these C stars are AGB stars undergoing 
the third dredge-up or, alternatively, that they underwent mass 
accretion from an old AGB companion, which subsequently 
evolved into a white dwarf. More exotic explanations will be 
discussed in Sect. 15.21 



2. Observations 

2.1. The target cluster: NGC 1978 

The globular cluster NGC 1978 is a memb er of the in termediate- 
age (t - I - 3 Gyr) cluster population ( Girardi et al . 1995) in 
the LMC. This population corresponds to one of the two peaks 
in the LMC metallicity distributi on at [Fe/Hj - - 0.31 (with a 
sprea d of tr = 0.15) reported b y lCole et al.l dSoOSl) . In a recent 
work. IMucciarelli et al.l (l2007al) indeed derived a metall icity of 
M/HJ - -0.37 that is based on the findings from Ferraro et al.l 
20061) . who claimed [Fe/H] ^ -0.38 + 02 an d that [a/H] is al- 
most solar In contrast to that. lHill et al.l (l2000l) quoted a grossly 
deviating value of [Fe /H] - -0.96 ± 0.20. They investigated two 
stars from this cluster (LE8 and LE9, see also Sects. [3]and|4]i and 
concluded that [O/Fe] = +0.37 ± 0.10. These find i ngs ha ve not 
been reproduced by other groups. IMucciarelli et al.l ( l2008l) rather 
found a mild depletion in alpha elements on average and as- 
cribe t he differences to Hill et al. (2000) to the deviating metal- 
licity. IMucciarelli et al.l (120081) . who performed detailed abun- 
dance studies of re d giants with UVES sp ectra, also did not con- 
firm the results of lAlcamo et al.l ( Il999l) . who concluded from 
multicolour CCD photometry that NGC 1978 consists of two 
sub-populations differing by 0.2 dex in metallicity. According to 
their work, the north-western half of the cluster is more metal- 
rich and younger than the south-eastern half (plus differences 
of the sub-populations in the helium content). This result has to 
be seen in conjunction with a peculiarity of the cluster, namely 
its high eUipticity of e = 0.3(0 which has been found inde- 
pe ndently by variou s group s and was most recently confirmed 
by IMucciarelli et al.l (l2007al) . A merger event that could explain 
both the elongated shape and the suspec ted two sub-popul ations 
was, however, ruled out for example by lFischer et al.l (119921) . A 



' The flattening of a globular cluster is defined as 6 = (a-b)/a, where 
a is the major axis and b is the minor axis of the cluster. 



2.2. Spectroscopy 

We recorded high-resolution near-infrared spectra of 12 AGB 
stars in the LMC globular cluster NG C 1978. Ten stars from our 
target list were already identified by ILlovd Evansl (1X980) as red 
giant stars (and are tagged LE in th is work) but largely without 
information about the spectral type. Froge l et al.l d 19901) listed a 
number of AGB stars in NGC 1978, comprising also the stars 
from LE, and gave information about the spectral type, the clus- 
ter membership, and near-infrared photometry data. Based on 
this information, we selected ten targets that could be observed 
with the P hoenix s pectrograph mounted at the Gemini South 
Telescope (iHinkle et al. 199 8, 2003) with reasonable exposure 
times. Additionally, we constructed a colour-ma gnitude diagram 
from Two Micron All Sky Survey (2MASS. 1skrutskie et all 
2006) data and picked two more stars that, judging from their 
K magnitude and (J - K) colour index, are also located on the 
AGB of NGC 1978 (stars A and B). Details about the observa- 
tion targets together with the observing log are given in Table [T] 
The distribution of all the targets within the cluster is illustrated 
in Fig.[T] 

Our observing programme required observations at two dif- 
ferent wavelengths, one in the H band and one in the K band. 
For this purpose, we utilised the Phoenix order sorting filters 
H6420 and K4220. However, poor weather conditions in the 
last of the 4 observing nights prohibited to record H-band spec- 
tra for the targets LE2 and LE7. For the star LE9, we ob- 
tained spectra from queue mode observations in the semester 
2008A. The exact wavelength set tings are similar to the ones 
described in lLebzelter et al.l (l2008l) . i. e. in the K band our spec- 
tra run approximately from 23 580 to 23 680 A. In the H band 
(15 530 - 15 600 A), the spectral region observed was shifted to 
slightly higher wavelengths to cover a larger part of the CO 3-0 
band head at 15 581.6 A = 6417.8 cm"'. The slit width was set 
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Table 1. NGC 1978 targets and log of observations (LE stands for lLlovd Evanslll980h 



ID 


Type 


J 


K 


RA 


Dec 


H band {A^ = 


15585 A) 


(S/N)h 


Kband(^, = 23670 A) 


(S/N)k 






2MASS 


J2000 












LE4 


M 


12.347 


11.199 


05 28 43.72 


-66 14 03.7 


2006 Dec 03: 


3 X 1000 s 


65 


2006 Dec 05: 3 x 1000 s 


65 


LE6 


C 


12.090 


10.707 


05 28 46.27 


-66 13 56.4 


2006 Dec 07: 


3 X 800 s 


70 


2006 Dec 06: 3 x 800 s 


55 


A 


M 


12.253 


11.116 


05 28 46.23 


-66 13 24.8 


2006 Dec 03: 


3 X 1000 s 


55 


2006 Dec 05: 3 x 1000 s 


65 


LE3 


C 


11.490 


9.676 


05 28 44.49 


-66 14 03.9 


2006 Dec 03: 


3 X 400 s 


105 


2006 Dec 05: 3 x 400 s 


120 


LEIO 


M 


12.917 


11.802 


05 28 44.44 


-66 13 59.9 


2006 Dec 03: 


3 X 1800 s 


60 


2006 Dec 05: 3 x 1800 s 


70 


LE5 


M 


12.527 


11.387 


05 28 43.64 


-66 13 53.0 


2006 Dec 03: 


3 X 1200 s 


50 


2006 Dec 05: 3 x 1200 s 


60 


B 


C 


12.367 


11.078 


05 28 43.65 


-66 14 09.6 


2006 Dec 07: 


3 X 700 s 


50 


2006 Dec 06: 3 x 1000 s 


55 


LEI 


C 


12.282 


10.467 


05 28 48.50 


-66 14 59.9 


2006 Dec 03: 3 x 600 s 


50 


2006 Dec 05: 3 x 600 s 


115 


LE2 


C 


12.967 


11.266 


05 28 48.62 


-66 15 18.7 








2006 Dec 06: 3 x 1200 s 


65 


LE7 


C 


12.764 


11.214 


05 28 47.85 


-66 14 44.0 








2006 Dec 06: 3 x 1000 s 


75 


LE8 


M 


13.018 


11.916 


05 28 48.47 


-66 14 38.7 


2006 Dec 03: 


3 X 2000 s 


30 


2006 Dec 05: 3 x 2000 s 


40 


LE9 


M 


13.315 


12.259 


05 28 50.73 


-66 14 44.2 


2008 Mar 19: 4 x 1800 s 


45 


2006 Dec 06: 3 x 2000 s 


50 



NGC 1978 

2MASS - K band 




N 

.J 

Fig. 1. Distribution of observed targets in the cluster NGC 1978. 
The picture was acquired from the 2MASS catalogue using 
Aladin ([Bonnarel et al. 2000). Numbers refer to the nomencla- 
ture of ILloyd Evansi (il980 ). The additional targets A and B 
were chosen according to their position in the 2MASS colour- 
magnitude diagram. The stars A, 4, 5, 8, 9, and 10 possess 
an oxygen-rich atmosphere while B, 1, 2, 3, 6, and 7 are car- 
bon stars. The cluster is elongated and its major axis stretches 
roughly from south-east to north-west. The displayed sky region 
is 5.'0 X 5'0 wide. The coordinates of the image centre (close to 
LE3) read a = 05^28'M5!00, 6 = -66°14'14'.'00 (J2000). 



to 0'.'35 (the widest slit) which resulted in a spectral resolution 
of/J = /l/A^ = 50000. 

The total integration time per target and wavelength setting 
ranged between 20 (LE3) and 120 (LE9) minutes. For each target 
we observed at two or three different positions along the slit. 
Each night, we also recorded a spectrum of a hot star without 
intrinsic lines in the respective wavelength regions in order to 
correct the spectra of our programme stars for telluric lines. The 
resulting signal-to-noise was from 50 to above 100 per resolution 
element (~ 4 pixels). See Table[T]for details. 



2.3. Data reduction 

The data reduction procedure was carr ied out as described at 
length for example in lSrnith et alj (|2002|) and in the Phoenix data 
reduction IRAF tutorial0 

To correct the AGB star spectra for telluric lines, we ac- 
quired the spectrum of an early type star without intrinsic lines. 
The telluric absorption features were removed using the IRAF 
task telluric in the K band. The H-band spectra are almost 
free of telluric lines, they were, however, also processed in 
the same way to remove the fringing. In the H band, we did 
the wavelength calibration for a K-type radi al velocity stan- 
dard (HD5457, Vr = 5.1kms''. rWUsonlll953h recorded along- 
side with the progr amme stars. Using the Arcturus atlas from 
iHinkle et al.l d 19951) . we identified several features (OH, Fe, Ti, 
Si lines) in the spectrum and derived a dispersion solution. The 
relation was then applied to the remaining spectra. This proce- 
dure allowed us to derive radial velocities from the H-band spec- 
tra. For the K-band spectra we did a direct calibration using the 
CO lines in the spectrum of an M-type target. That solution was 
then also applied to the carbon-star spectra. 

3. Data analysis 

3. 1 . Contents of the observed wavelength ranges 

The wavelength ranges were chosen such that from the H band 
spectra we could derive the stellar parameters and the C/O ra- 
tio. Subsequently, we inferred the carbon isotopic ratio '^C/'^C 
from the K band. In practice, the parameters are tuned iteratively. 
The region in the H band is widely used i n the literature (e .g. 
ISmith et al.ll200l iMcSavenev et~al]l2007t lYong et al.ll2008h to 
derive oxygen abundances from the contained OH lines. From 
the relative change of those features in comparison to the band 
head from the '^C"'0 3-0 vibration transition, it is possible to 
determine the C/O ratio in the stellar atmosphere. A weak CN 
line and a few metal lines (Fe, Ti, Si) aid in the pinning down of 
the parameters, especially the effective temperature. 

The K-band spectra comprise a number of CO lines from first 
overtone (Av - 2) transitions. Beside features from the main iso- 
topomer '-C'^O, we also find some '^C'^O lines in this region. 
We derived the carbon isotopic ratio by fitting these lines. Apart 

- Available at http://www.noao.edu/usgp/phoenix/ IRAF is 
distributed by the National Optical Astronomy Observatories, which 
are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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from the CO lines, there is also a single HF line blended with a 
'^CO feature. 

In Fig. |2] we give an overview about our observations. For 
the oxygen-rich stars, the key features in the H band are marked. 
In the K band, we indicate the position of the ' ^'CO lines. The se- 
lection of the wavelength ranges was driven by the oxygen-rich 
case. From Fig.|2]it is obvious that the regions are not well suited 
for the analysis of carbon-star spectra. Both H- and K-band spec- 
tra are crowded with features from the CN and C2 molecules, 
occurring in addition to the CO lines. The polyatomic molecules 
C2H2, HCN, and also C3 contribute to the opacity by means of 
many weak lines that form a pseudo-continuum (see Sect. |3.3.2] i. 
The major part is due to C2H2 absorption that increases with de- 
creasing effective temperature. In the H band, the CO band head 
is strongly affected by neighbouring features. Generally speak- 
ing, there is not a single unblended line in the carbon-rich case. 
The situation in the K band, however, is not that bad. The '^CO 
lines are still visible, although blended with CN features, but 
there are only a few C2 lines. The pseudo-continuum is essen- 
tially all due to C2H2, contributions from other molecules are 
negUgible. 

3.2. Synthetic spectra 

We compare our observations with synthetic spectra based on 
model atmospheres that were calculated with the COMARCS 
code (details will be given in a forthcoming paper by 
lAringer et af] |2009|) . COMARCS is a modified version 

of the MARCS co de (originating from iGustafsson et alj 

1975; J0rgensen et al. 1992, recently extensively described by 
Gustafsson et al. ,2008) . In the model calculations, we deduced 
the temperature and, accordingly, pressure stratification assum- 
ing a spherical configuration in hydrostatic and local thermody- 
namic equilibrium (LTE). For stars showing only a small vari- 
ability, LTE is a reasonable approximation. Possible deviations 
from LTE have a smaller influence on the abundance determi- 
nation than larg e deviations fr om the hydrostatic equilibrium 
(see for example IjohnsonlFl 9941 for a discussion of non-LTE ef- 
fects in cool star atmospheres). In the evaluation of the chemical 
equilibrium, which is a consequence of LTE, we take all rele- 
vant opacity sources into account in both model calculation and 
spectral synthesis. In the oxygen-rich case H2O, TiO, CO, and 
CN are major opacity contributors, while HCN, C2H2, C2 (and 
others) are important in the carbon-rich case. The opacity coef- 
ficien ts utilised by COMARCS are calcu lated with the COMA 
code (lAringeij|2d00l; |Aringer et al."2009'). The atomic line data 
were taken from VALD (Kupka et al. 2000), an overview of the 
molecular l ine lists used wi th all the sources documented can 
be found in lLederer & Aringer (2009). Model atmospheres and 
synthetic spectra that were calculated following the method out- 
lined above have been sh own to de scribe the spectra of cool giant 
stars appropriately (e. g. iLoidl et al. 2001; Aringer et al. 20021). 
We applied the same procedure successfully in paper I for our 
analysis of the AGB stars in the LMC cluster NGC 1846. 

From the parameters that determine an atmospheric model, 
we held the mass and the metallicity constant. The respec- 
tive values are M - 1.5 Mq and [M/H] - -0.4, and 
were taken from the literature (see Introduction). In the model 
calculations, the microturbulent velocity was set to ^ - 
2.5 km s ', which is in t he range that is found for atmospheres 
of low-mass giants (e. g. Aringe r et al.ll2002t ISmith et al.ll2002l: 
[Gautschy-Loidl et al. 2004). By varying the remaining parame- 
ters effective temperature T^ff [K], logarithm of the surface grav- 
ity log(§[cms"^]), and carbon-to-oxygen ratio (C/O), we con- 



structed a grid of model atmospheres and spectra. This was done 
such that we cover the Tetf and log g range resulting from colour- 
temperature relations and bolometric corrections applied to our 
sample stars. The step size in effective temperature was set to 
50 K, whilst the surface gravity was altered in steps of 0.25 on 
a logarithmic scale, i. e. for \ogg, ranging from 0.0 to H-0.5. For 
some carbon stars in our sample we got \ogg < 0.0 from the 
colour transformations and the adopted mass. The spectral fea- 
tures used in the analysis show only a minor dependence on this 
parameter thus we fixed \ogg = 0.0 in the analysis of the carbon 
stars. In this way we also avoid convergence problems for the 

atmospheric models. 

The element composition is scaled solar jGrevesse & Noelsl 

mH), but we assumed an oxygen over-abundance of -1-0.2 dex 
motivated by the results of our work on NGC 1 846 (paper I) and 
by the paper of lHill etall (|2000). We altered the C/O ratio in 
steps of 0.05 in the oxygen-rich case and 0.10 in the carbon-rich 
case. This was done by changing the carbon abundance while 
leaving the other abundances untouched. 

The synthetic spectra cover a wavelength range of 6400 - 
6450cm-' (15 504 - 15 625 A) in the H band and of 4215 - 
4255 cm^' (23 500 - 23 725 A) in the K band. The spectra were 
first calculated with a resolution of R - A/AA = 300000 and 
then convolved with a GaussiarQ to match the resolution of our 
observed spectra (R = 50000). By applying another convolution 
with a Gaussian profile we took the macroturbulent velocity into 
account. To determine the carbon is otopic ratio '^C/'^C (whic h 
is about 89.9 in the Sun according to lAnders & Grevessdl 19891) . 
we altered this parameter as well. The carbon isotopic ratio has 
virtually no effect on the model structure and was consequently 
only considered in the spectral synthesis calculations. 

3.3. Determination of abundance ratios 

We derive the effective temperature T^ff, the surface gravity 
log g, the C/O and '^C/'^C ratio of our targets by fitting synthetic 
spectra to the observations. The initial guesses for Teff and log^ 
were obtained using colour-temperature relations and bolomet- 
ric corrections from the literature. We converted the /r2MASS and 
(J - g) 2MASS values into the Johnson system with the formulae 
given m ICarpented(|200 l'). For the oxygen-rich targets, we esti- 
mated T^g and log g from the relations given by Houdashel t et al.l 
(200^ In the case of the carbon stars, we utilised the work of 
Wood et all (Il983h and lBesseU et all (11983 ) to obtain bolomet- 
ric corrections (BC^:) and effective temperatures, respectively. 
From the K magnitude together with BCa: and the distance mod- 
ulus (see Sect. |2]i we derived log L/Lq, and subsequently log^ 
with the obtained Tetf and by assuming M = 1.5 Mq. The 
reddening of this cluster was take n into account, but it is low 
(E(B - y) = 0.05 and Ak = 0.017. iGrocholski et al.ll2007l) and 
has no significant influence on the transformations. 

For each target we fit the parameters in a two-step process. 
The idea is to fit the spectral ranges at once rather than tuning 
individual abundances (except for carbon) to fit single spectral 
features. From the H band, we could in this way constrain the 
stellar parameters (Teff and \ogg) and the C/O ratio. Variations 
in log g had the smallest effect on the synthetic spectra. A small 
uncertainty in the stellar mass or radius estimate has thus only 
a minor influence on the derived abundance ratios. Using the 

^ The Gaussian has the functional form G(AA,cr) = 
exp{-(A/l/cr)V2} /(cr V27r), whereby cr^ = A/(2R) to reduce the 
effective resolution to a given R, or o",,, = (/l/2)(vt/c) to account for the 
macroturbulent velocity Vf. 
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LE5 



LE4 



LE6 



LE7 



LE2 



LE3 
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"CO '^CO '^CO+HF 
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'CO 2-0 P5 



'CO 2-0 P6 
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'CO 2-0 P8 
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Wavelength [A] 



15590 



15600 23580 23590 23600 23610 23620 23630 23640 23650 23660 23670 
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Fig. 2. Overview about observations: H-band spectra are shown on the left side, K-band spectra on the right. The spectra are ordered 
by increasing infrared colour A')2mass (from top to bottom), respectively. Some features are identified for the oxygen-rich targets 
(top), for the carbon-star spectra we only indicate the position of the low-excitation CO lines. Data of the target LE8 are not shown 
as they are of poor quality and we did not derive abundances from them. H-band spectra for the targets LE2 and LE7 could not be 
taken due to poor weather conditions. See text for details. 



parameters as obtained from the H band, we then calculated K- 
band spectra with a varying '^C/'^'C ratio to determine this pa- 
rameter as well. The carbon isotopic ratio does not influence the 
appearance of the H band spectrum significantly, as was verified 
by test calculations. 

To match the shape of the spectral features we also had to 
assume a macroturbulent velocity which reduces the effective 
resolution of the spectra. This artificial broadening includes the 
instrumental profile, but it is also used to imitate at least partly 
the dynamical effects in the stellar atmosphere, which become 
increasingly important for carbon stars. This is also why the 
adopted values were generally lower for the oxygen-rich tar- 
gets: they were fitted by applying a macroturbulent velocity of 
Vt = 3kms"'. The carbon stars displayed a stronger broadening 
of the features (cf. Fig.|2]l. We needed values of Vt - 8kms"' 
and even Vt = lOkms"' (LE6) to fit the spectra (compare 
also Ide Lavernv et al.ll2006l who find comparable values in their 
study). We consistently applied the same value for the H and the 
Kband. 

We tried to objectify the search for the best fit by applying a 
least-square fitting method. Although we were in this way able 
to narrow down the possible solutions to a few candidate spec- 
tra quickly, the final decision about our best fit was done by vi- 



sual inspection. The reason is that spectral features with a false 
strength or position (both due to imperfect line data) in the syn- 
thetic spectrum, or spectral regions with a high noise level dom- 
inate and confuse an automatic minimisation algorithm. As 
the last step in the fitting procedure is done by eye, we assign a 
formal error to the derived parameters. In the case that the best 
fit parameters lie in-between our grid values, we quote the arith- 
metic mean of the candidate models as our fit result. The error 
bars were estimated from the range of parameters of the model 
spectra that still gave an acceptable fit. 

The formal errors for the derived C/O and '^C/'^C ratios 
given in Table [3] compare well with the uncertainties found by 
some basic error estimates, as will be shown in the following. 
Therefore the quintessence of the discussion in Sect. |5] is on a 
sound footing. We consider typical uncertainties in the stellar 
parameters and investigate the influence on the derived abun- 
dance ratios. The correlations between the stellar parameters and 
possible systematic errors from the model syntheses have not 
been taken into account in the error analysis. For the oxygen- 
rich case, we start from a baseline model with T^g - 3750 K, 
Vt = 3kms \ logg = 0.25, and C/O = 0.15. Changes of 
AFeff - 100 K, Alogg = 0.25, and Avt = 2kms~' correspond 
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to changes in A(C/0) of 0.02, 0.03, and 0.01, respectivel>Q. 
Summed in quadrature this results in a typical total uncertainty 
of 0.04 for the C/O ratio of our M-type targets. The same ex- 
ercise for the carbon isotopic ratio results in A('^C/'''C) of 2, 
2, and 3 (parameter dependence as above), and thus a total un- 
certainty of 4. Concerning the carbon-rich case, we start from a 
carbon-rich model with parameters T^ff = 3350 K, \ogg - 0.25, 
Vt = 8kms"', C/O = 1.40, ^^C/"C^115 and vary the pa- 
rameters as described above. This leads to a total A(C/0) - 
VO.l^ + 0.05^ -1- O.P = 0.15. For the carbon isotopic ratio we 
find AC'^C/'^C) = VTPTTs^TsO^ ^ 37. 

A considerable part of our discussion will be concerned with 
the '^C/'^C ratios of our targets. Since the derived isotopic ra- 
tios d epend on the parameters of the '^CO lines (taken from the 
list of iGoorvitoh & Cha ckerian 1994, see also AppendixlAl. we 
want to assess possible systematic errors in the line strengths. In 
the case that the predicted line strengths in the list are too large, 
one would have to increase '^C/'^C (equivalent to a decrease 
of the '^C abundance) in order to fit a '^CO feature compared 
with the case where the predicted strengths are correct. The ac- 
tual '^C/'^C ratio would then be lower than the derived value. 
For some oxygen-rich stars in our sample we get carbon isotopic 
ratios that are close to the value of the CN cycle equilibrium (4- 
5). A necessary further reduction is thus not very likely. If the 
theoretical line strengths are too low, the above argument is re- 
versed. One would assume a lower '^C/'^C ratio to fit the '^CO 
features, which means that one would underestimate the actual 
isotopic ratio. Since we, of course, use the same line list for the 
analysis of both M and C stars, this would shift up the '^C/'^C 
for all targets in Figs.|2l|8] and|9] Similar to a \ogg uncertainty, 
the values for the carbon stars would be stronger affected. Hence, 
underestimated '^CO line strengths would relax the necessity of 
efficient additional mixing processes on the RGB as discussed in 
Sect.ia 

In the next section, we discuss M and C stars separately. The 
features contained in the spectra depend on the chemistry regime 
and some issues cannot be discussed in a general way. 

3.3.1. M-type stars 

In the search for a good fit we took advantage of the way features 
react to parameter changes. An increasing temperature weakens 
all features in the H band. The OH/Fe blend at approximately 
15 570 A is a good temperature indicator The two neighbour- 
ing OH lines and the CO band head react less on temperature 
changes. The remaining lines show only a weak dependence on 
Jeff- The K band is largely insensitive to temperature changes, 
only the '^CO lines show a weak dependence on temperature 
which adds to the uncertainty in determining the isotopic ra- 
tio. The only feature strongly reacting when altering T^ff in the 
model calculations is the blend including the HF line, its strength 
decreases when the temperature is increased. 

An increase in C/O enhances the strength of the CO band 
head and causes stronger CN lines, while the OH lines get 
slightly weaker, especially the feature blended with an Fe line. 
Lowering C/O decreases the strength of the CO features in the 
K band. 

The CO band head is also sensitive to changes in \ogg, 
while the other features practically are not. Changes in the sur- 

Note that we also include changes of the macroturbulent velocity 
in the error analysis. The equivalent width of a spectral line is unaf- 
fected by Vt, but in the spectrum synthesis method uncertainties in this 
parameter influence the derived abundances. 



face gravity also affect the measured carbon isotopic ratios (see 
Sect. |3.3l l. Two of the '^CO lines in the spectra of the oxygen- 
rich stars are unblended ('^CO 2-0 R18 and '^CO 2-0 R19). 
Blended '^CO features were used to check the isotopic ratio de- 
rived from the clean lines. 

The continuum is well defined in the oxygen-rich case. We 
utilised the least-square method in the process of finding the best 
fit. 

3.3.2. C-type stars 

Several things make the fitting of carbon stars more difficult 
compared with the case of M-type stars. Foremost, the quality 
of the line data for the molecules appearing in the carbon-rich 
case hampers the qualitative analysis profoundly. For CN, C3, 
and C2H2 we used the SCAN database (J0rgensen 1997). The 
line positions in the computed CN list are not accurate enough 
for modelling high-resolution spectra. With the h elp of mea- 
sured line positions for CN that were compiled bv iDavis et al.l 
(l2005h . we were able to correct the wavelengths of the strongest 
lines. While the results were satisfying in the oxygen-rich case, 
in the carbon-star spectra many additional weak lines show up 
that coul d not be corrected . The case of C2 (we used the line 
data from lOuerci et al.lI974l) . producing a wealth of spectral fea- 
tures, is even more problematic: no observed reference line list is 
available for this molecule, so both line positions and strengths 
are subject to large uncertainties. We corrected the line list man- 
ually by removing strong features that did not appear in any of 
our observations. Several features were shifted to other wave- 
lengths where it was evident from the observations that the lines 
are at the wrong position. Unlike for M stars, the H-band spec- 
tra of carbon stars are also affected by the carbon isotopic ra- 
tio. Lines from '^C'^C or '^C''*N are important in some blends, 
however, the quality of the line data for these features could not 
be assessed for the above described reasons. 

The molecules C2H2 and C3 are incorporated into our cal- 
culations via low-resolution opacity sampling data. The absorp- 
tion thus becomes manifest as a pseudo-continuum in the spec- 
tra. A possible occuiTence of strong lines or regions with par- 
ticularly low absorption cannot be reproduced with this ap- 
proach. The pseudo-continuum level reacts sensitively to tem- 
perature changes. In general, an increase in Teff reduces the fea- 
ture strength. The relative changes of different line strength can 
be used to constrain the temperature range. The lower the tem- 
perature gets, the higher the contribution of C2H2 and C3 be- 
comes, whereas C2H2 dominates the absorption. An increase in 
the C/O ratio has the same effect. The implications for abun- 
dance determination of the ill-defined continuum in carbon stars 
is discussed in detail in paper I. We want to stress here the 
consequences for the method. The absorption in the pseudo- 
continuum varies with temperature requiring scaling of the ob- 
served spectrum before a comparison. Due to this scaling the 
value of changes in the same manner, so that model spectra 
with lower effective temperatures (causing a lower flux level due 
to an increased pseudo-continuum) always result in a lower 
pretending to fit the observations better. This is of course an arti- 
ficial effect, and thus we cannot rely on the least-square method 
as objective criterion. 

A change in log g has only small effects on the carbon-star 
spectra. A higher surface gravity produces a higher pseudo- 
continuum in the K band. In general log g only marginally affects 
the spectral lines. However, in the determination of '^C/'^C the 
uncertainty in log g has to be taken into account. In the K band, 
all '^CO lines are blended with other features. The carbon iso- 
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topic ratio is already so high that the strengths of the CO lines 
have become rather insensitive to changes in this parameter As 
a consequence, even the small uncertainties in log g correspond 
to large changes in the isotopic ratio, which adds to the errors. 

Variations of the C/O ratio have the strongest impact on the 
H-band spectra when C/O is slightly above 1. The strength of 
the CO band head rapidly drops when C/O is increased to 1.3, 
approximately. Then a saturation sets in and the strength of the 
CO band head varies slowly with C/O. This behaviour and strong 
CN features sitting in the band head hamper an accurate C/O 
determination. The lines of C2 and CN increase in strength when 
C/O rises, both in the H and the K band. The CO lines in the K 
band decrease in strength for higher values of C/O. The changes 
are, however, small and do not allow for a determination of the 
C/O ratio from the K-band spectra alone. 

Figure[3]is, as an example, a fit for the star B, which is a car- 
bon star. The star has an effective temperature of about 3350 K 
(see also Table [3] for the other fit parameters). The pseudo- 
continuum contribution is relatively weak, although there are 
almost no line-free regions in both spectral bands. A number 
of features is successfully reproduced by our models, the devi- 
ations in other regions are most probably due to uncertain line 
data. 

To measure the '^C/'^C ratio in the carbon-star spectra, we 
utilised the '^^CO lines at 23 579.9, 23 590.2, and 23 601.0 A (the 
leftmost lines indicated in the right panel of Fig.|2|. In the vicin- 
ity of the other '"'CO lines, there are obviously absorption fea- 
tures missing in our synthetic spectra (see Fig. [3]). We cannot rule 
out that our analysis of a limited number of line blends could in- 
troduce a systematic eiTor in the infeiTed carbon isotopic ratios. 
Therefore, we emphasise that it would be worthwhile to reassess 
the isotopic ratios in other wavelength regions, particulary in the 
light of the discussion in Sect. |5] 

We also identified four low-excita tion '^CO 2-0 lines in the 
K band (see also Fig.|2|. According to lHinkle et alj (|1982|) . low- 
excitation lines form in the outer atmospheric layers, opposite 
to high-excitation lines that show characteristics similar to sec- 
ond overtone transitions (Ay = 3). In all the carbon-star spectra, 
these lines could not be fitted with our synthetic spectra. The 
calculations underestimate the line strength, suggesting that the 
extended atmospheres of luminous carbon stars are not well re- 
produced by our hydrostatic models. 



3.4. Stellar evolutionary models 

The stellar evolutionary models presented in this paper hav e 
been computed with the FRANEC code jChieffi et al.1 Il998h . 
The release we are currently using is optimised to properly com- 
pute low-mass models during their AGB phase. Up-to-date in- 
put physics, such as low-temperature carbon-enhanced opaci- 
ties, have been adopted ( Cristallo et al. 2007; Lederer & Aringer 
12009 ). Physical phenomena, such as hydrodynamical instabil- 
ities at radiative-convective interfaces a nd the mass-loss rate, 
have been prop erly taken into account (IStraniero et al.l [20061; 
ICristallo et al.l 1200 9). The inclusion of an additional mixing 
mechanism taking place below the convective envelope (usually 
referred to as extra-mixing) during the RGB and the AGB phase 
was described in paper I. 



Table 2. Radial velocities of NGC 1978 targets 



ID 


vv[kms '] 




A 


+295.6 


t 1.0 


LE4 


+296.4 


t 0.8 


LE5 


+292.9 


t 0.9 


LE8 


+290.4 


t 1.3 


LEIO 


+292.5 : 


t 0.7 


B 


+290.6 


t 1.3 


LEI 


+297.5 


t 1.6 


LE3 


+291.5 


t 1.9 


LE6 


+290.2 


t 1.7 


mean 


+292.9 


t 1.3 


velocity dispersion (tr) 




2.6 



4. Results 

4.1. Cluster membership 

We derived the heliocentric radial velocity by cross-coiTelating 
the H-band spectra with a template spectrum (using the IRAF 
task fxcor). Data for the star LE9 were taken in the semester 
2008A, and no radial velocity standard (K-type star) was 
recorded, so no information could be deduced for this target. 
While the quality of our LE8 data is too low to derive abundance 
ratios, it is still adequate to measure the radial velocity (only 
with a slightly larger eiTor than for the other targets). The sam- 
ple shows a relatively homogeneous radial velocity distribution, 
the spread is quite narrow, so we conclude that all our targets are 
actually cluster members. The results are summarised in Table|2] 
The value Vr = +292.9+ 1 .3 km s"' (with a velocity dispersion of 
cr,, = 2.6 km s ') we find is well in line w ith earlier determina- 
tions (see Introduction). As a comparison. ICarrera et a n (I2008h 
measured radial velocities in four fields with different distances 
to the LMC centre. The mean values found in the individual 
fields range from +278 to +293 km s"'. The velocity dispersion 
in the metallicity bin with [Fe/H] > -0.5 is o",, = 20.5 km s"'. 

The errors for the carbon stars are systematically larger than 
for the oxygen-rich stars. This can easily be understood by look- 
ing at Fig.|2] The features in the C-rich case are usually broader 
and consequently the peak in the cross-coiTelation function (and 
thus the FWHM) is broader, too. 

4.2. Stellar parameters and abundance ratios 

We summarise our fit results in Table [3] In the first column, 
we list the star identifier (compare Fig. [U, followed by the 
K magnitude and the colour index (J - K) taken from the 
2MASS database. In the next two columns, we quote the effec- 
tive temperature and luminosity (rounded to 100 Lq) resulting 
from the respective colour calibrations and bolometric correc- 
tions (see Sect. I3.3l l. The parameters resulting from our fitting 
procedure (Teff, \ogg, Vt, C/O, '^C/'^C) are listed in the subse- 
quent columns. We note that none of our target stars is a large 
amplitude variable, thus the influence of variability on the stel- 
lar parameters can be safely ignored. A detailed discussion of 
the small amplitude variability (1.5 mag in R) will be given else- 
where (Wood & Lebzelter, in preparation). 

The derived effective temperatures are systematically higher 
than the T^s values inferred from the infrared colour transfor- 
mations (compare the resulting Teff^c to T^^i deduced from the 
fit in Table |3]l. Better agreement between the two temperature 
scales is found when we adopt a scaled solar oxygen abundance 
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Wavelength [A] 

Fig. 3. Observation and fit of the carbon star B in NGC 1978. The fit parameters are given in Table |3] The relatively high effective 
temperature and low C/O ratio compared with the other carbon stars in our sample allowed a reasonable fit. Deviations can to a high 
fraction be ascribed to uncertain line data. The strong low-excitation CO lines in the K band are not reproduced by the models. See 
text for details. 



Table 3. Data and fit results for oxygen-rich (first group) and carbon-rich (second group) targets in NGC 1978 



ID 


K 


J-K 










V, 


C/O 


'^c/'-*c 




2MASS 


[K] 




[K] 




[kms-'] 






A 


11.116 


1.137 


3600 


5300 


3825 


0.50 


3 


0.23 ± 0.05 


16 ± 3 


LE4 


11.199 


1.148 


3600 


5000 


3725 


0.38 


3 


0.18 ±0.03 


13 ±4 


LE5 


11.387 


1.140 


3600 


4200 


3775 


0.25 


3 


0.18 ±0.05 


6±2 


LE9 


12.259 


1.056 


3800 


2000 


3900 


0.38 


3 


0.13 ±0.03 


12 ±2 


LEIO 


11.802 


1.115 


3700 


3300 


3900 


0.38 


3 


0.18 ±0.05 


9±3 


B 


11.078 


1.289 


3176 


4400 


3350 


0.00 


8 


1.35 ±0.10 


175 ± 25 


LEI 


10.467 


1.815 


2554 


6200 


2600 


0.00 


8 


> 1.50 


> 200 


LE3 


9.676 


1.814 


2556 


12900 


2600 


0.00 


8 


> 1.50 


> 200 


LE6 


10.707 


1.383 


3043 


5800 


3100 


0.00 


10 


1.30 ±0.10 


150 ± 50 



rather than an over-abundance of 0.2 dex, which was our stan- 
dard choice in the analysis (we refer to Sects. [1] and 13.21 for 
details). However, this is n ot too surprising s ince t he colour- 
temperature relations from iHoudashelt et al.l (l2000l) were de- 
rived from synthetic spectra based on s caled solar abunda nces. 
We want to mention that the results from Smith et al. (2002, their 
Fig. 8) imp ly [O/Fe] < 0.0 fo r the LM C, contrary to our assump- 
tion. Also, iMucciarelli et al.l (l2008l) found that the other alpha 
element abundances are roughly scaled solar ( [g/Fe] ^ 0.0) in 
NGC 1978. However, in another recent work i Goudfrooii et al.l 
(|2009) derived [a/Fe] = 0.2 for the LMC cluster NGC 1846. 
Apart from the influence on the temperature scale, a higher oxy- 
gen abundance would shift the derived C/O ratios only to slightly 
higher values (Fig. |4]i. The derived carbon isotopic ratios are 
marginally affected by an oxygen over-abundance. 

We found that for the five oxygen-rich stars within our sam- 
ple the C/O ratio is varying little, the values range from 0.13 to 
0.18 with a typical uncertainty of +0.05. The isotopic ratios vary 
in the range between 9 and 16 with an uncertainty ranging up to 
±4. Considering the error bars, all M stars occupy more or less 
the same region in Fig. |5] where we display the measured C/O 
and '-C/'^C of our sample stars. No target is offset from the spot 
marking the abundance ratios expected after the evolution on the 



first giant branch (C/O - 0.2, ^^C/^^C- 10). In other words, 
we do not find conclusive signs of third dredge-up among the 
oxygen-rich stars in the cluster. This is also consistent with the 
results from Llovd EvansI (Il983h who did not find S-type stars 
in NGC 1978. The total number of our targets is small, so gen- 
eral statements based on our results are rather weak due to the 
poor statistics. However, the fact that we find oxygen-rich stars, 
in which the TDU does not seem to be active, does not exclude 
that this phenomenon is at work in other stars of the cluster 

In fact, we also identify a sub-sample of carbon-rich stars 
belonging to NGC 1978. Satisfactory fits could only be achieved 
for the two hottest carbon stars (B and LE6) in our sample. We 
also derived C/O and '^C/'^'C. The error bars are larger than 
for the M stars, the reasons for that are outlined in the previ- 
ous sections. For the two cool objects LEI and LE3, we could 
only derive lower limits for C/O and '^C/'^C. For an increas- 
ing carbon content, the features saturate, thus we cannot give a 
reliable upper limit for C/O. The C/O ratios th at we found are 
in line with the results obtained by other groups. iMatsuura et al.l 
(2005) adopted C/O > 1.4 to explain molecular features in low- 
resolution spectra of LMC carbon stars. Investigations of plan- 
etary nebulae in the LMC exhibit a range of C/O ratios from 
slightly above 1 up to even 5 (e. g. IStanghellini et a"Lll2005l) . We 
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Observation: NGC 1978 - LE5 
Model: [O/H]=0.0, Tg|,=3700 K, log g=0.25, C/O=0.20 
Model: [O/H]=0.2, Tg|f=3800 K, log g=0.25, C/O=0.15 
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Fig. 4. Fit of the H band spectrum of the target LE5 adopting different oxygen abundances. A value of [O/Fe] = 0.2 (which we used 
in our analysis) leads to a higher temperature and a lower C/O ratio compared with a scaled solar oxygen abundance ([O/Fe] = 0.0). 
The two model spectra both give reasonable fits to the observations, thus a decision about the actual oxygen abundance cannot be 
made from our data. However, lower effective temperatures (and thus the "no oxygen over-abundance" scenario) agree better with 
the values derived from colour-temperature relations. 
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Fig. 5. Measured abundances ratios C/O and '^C/^^C for the 
targets in NGC 1978. Oxygen-rich targets appear to the left of 
the dotted line at C/O = 1, carbon-rich targets are on the right. 
The box in the upper right corner marks the lower limits for the 
targets LEI and LE3. There are no obvious signs of TDU among 
the M stars (see left panel of Fig.|2]for a zoom to the M-star data 
points), S stars are completely lacking in our sample of cluster 
stars. Additionally, from the carbon-star data we cannot detect a 
saturation level for the carbon isotopic ratio which would hint at 
the occurrence of extra-mixing. 



are not aware of measurements of the carbon isotopic ratio di- 
rectly from carbon-rich AGB stars in the LMC o ther than the 
ones presented here and in paper I. iRevniers et al.l (2007) anal- 
ysed the post-AGB object MACHO 47.2496.8 ([Fe/H] = -1.4) 
in the LMC and derived C/O > 2 and '^c/'^c ~ 200. A pe- 



culiar combination of C/O and '^C/'^C values even more ex- 
treme than those of our carbon-rich targets was presented by 
Ide Lavernv et all (12006). They found C/O < 1.2 and ^^C/"C > 
300 for BMB-B 30 (with [M/H] = -1.0) in the SMC. 

Qualitative estimates based on the K-band spectra of the tar- 
gets LE2 and LE7 suggest that those have an effective temper- 
ature between the groups LE1-LE3 and B-LE6 (see Fig. |2]i. 
Except for the target B, the difference between the derived ef- 
fective temperature and the T^ff value obtained from the colour 
transformation is not quite as high as in the oxygen-rich case 
(about 50 K). We repeated the analysis for the star B with a 
scaled solar oxygen abundance (as for LE5, Fig. |4]i and also 
found a lower Teff, whereas the C/O ratio remained basically 
unchanged. Almost all spectral features in the carbon-rich case 
react sensitively to temperature changes, therefore Teff can be 
better defined than the C/O ratio. A remarkable result is that — 
opposite to our findings for NGC 1846 (paper I) — we do not find 
a saturation level for the carbon isotopic ratio. This indicates that 
extra-mixing is not at work in NGC 1978. 

In Fig. |6] we show a CMD of NGC 1978 (see Sect. |5] for a 
discussion of the photometric errors). We do not plot the bolo- 
metric magnitudes on the ordinate because of the inconsistent 
bolometric corrections applied to the targets. The CMD contains 
all targets that lie within a radius of 4.'0 with respect to the clus- 
ter centre (cf. Fig. [Ij. Data points with superimposed symbols 
refer to stars contained in our sample. The oxygen-rich stars 
(empty squares) seem to form a sequence with a C/O ratio in- 
creasing with luminosity which might hint at a mild dredge-up. 
However, this could be misleading, as we have to take the er- 
ror bars quoted in Table |3] into account. Moreover, the isotopic 
ratios do not appear in an ordered sequence. The carbon stars 
(red triangles) are brighter in K and have significantly redder 
colours. For the stars LEI and LE3, having the highest (J - K) 
values, the indicated numbers refer to lower limits for the abun- 
dance ratios. The two carbon stars with the lowest K magnitude 
are LE2 and LE7 which were excluded from our analysis. Both 
have luminosities comparable to the brightest oxygen-rich stars. 
A possib le explanation is tha t these stars are in a post-flash dip 
phase (ilben & Renzinil[T983h . displacing them about 1 magni- 
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Fig. 6. Colour-magnitude diagram based on 2MASS data for 
NGC 1978 (complete within a radius of 4.'0 with respect to the 
cluster centre). Empty squares refer to oxygen-rich objects in 
our sample, filled triangles depict carbon stars. Compare also 
Table |3] Details are discussed in the text. 



tude down from the average luminosity in this evolutionary stage 
in the CMD (see also the discussion in paper I). 



5. Discussion 

We seek to explain the following observed features of 
NGC 1978: 

1. The cluster harbours carbon stars. 

2. The M stars do not show conclusive signs of third dredge-up. 

3. There is no saturation of the carbon isotopic ratio. 

Except for the first point, this is more or less the opposite of 
what we found for the cluster NGC 1846 (paper I), the composi- 
tion of which is similar to the one of NGC 1978, even if the GC 
is slightly younger and, as a consequence, harbours more mas- 
sive AGB stars (with a mass of about 1.8 Mq). The M stars in 
NGC 1846 span a range in C/0 from 0.2 to 0.65, and ^^C/^^C— 
initially rising together with C/0 — does not exceed 60 to 70, 
even in the carbon stars. 

The first point in the list is a clear indication of ongoing TDU 
in the cluster The second point may at first seem to contrast this 
statement, while low number statistics could resolve this appar- 
ent contradiction; the lifetime of a thermally pulsing AGB star 
is short, so in a sample of a few stars we do not necessarily have 
to find a star with an enhanced carbon abundance. This possibly 
explains the large gap between the M and the C stars in terms of 
C/O (see Fig.|5]). An alternative scenario where the TDU is so ef- 
ficient that the star becomes carbon-rich after a sin gle dredge-up 
episo de occurs only at much lower metallicity (e. g. lHerwig et alj 
12000 ). 

Let us put aside the C/O error bars for a moment and spec- 
ulate about the group of M stars. The 2MASS data come with 
uncertainties in K of lower than 0.05, except for LEIO where 



the error is 0.0760 In any case, the luminosity sequence that 
the stars form is conserved. The star LE9 is the faintest at K, 
with iuk - 12.3, in the samp le and also the star with the lowest 
C/O ratio. ICioni etal.l(l200d) identify the RGB tip in the LMC 
at K^) - 12.1 mag. According to this value, LE9 could as well be 
an RGB star or located at the early- AGB. The targets LE5 and 
LEIO are possibly at the very beginning of the thermally pulsing 
asymptotic giant branch (TP- AGB) star phase. The case of LE4 
and, in particular, A is more difficult. Comparing these stars with 
those of similar luminosity in the sample of NGC 1846, we note 
that the latter are classified as S stars and present C/O ratios of 
0.4 or even higher. The C/O ratio we find for the star A (0.23) is 
slightly higher than the average of the other stars, perhaps com- 
patible with just one (the first) TDU episode. The bolometric 
magnitude of A also roughly corresponds to the limit we found 
for the onset of the TDU in NGC 1846 (mboi = 14.1 mag). 

5.1. Stellar evolutionary models 

In this section we present theoretical tracks computed with the 
FRANEC code, and we compare them with observations. In the 
left panel of Fig. |7] we report the evolution of the '^C/'^C ver- 
sus C/O curves from the pre-main sequence up to the early- AGB 
phase in the region where the M stars of NGC 1978 lie. In prac- 
tice, the end point of each curve (marked by a symbol in the left 
panel of Fig. |7]) represents the value attained at the tip of the 
RGB phase, which is conserved up to the onset of the first ther- 
mal pulse. In the right panel, instead, we extend the axes to also 
include the C-rich stars. We plot the whole AGB evolutionary 
tracks and mark the values attained after each TDU episode. 

We firstly concentrate on the O-rich stars of our sample. 
We carry on our analysis under the assumption that these stars 
have not experienced TDU, because they still are at the begin- 
ning of the TP- AGB phase or because they arrived on the AGB 
with a too small envelope mass. In this case, their initial surface 
composition has been modified by the occurrence of the First 
Dredge-up (FDU) only and, eventually, by an additional slow 
mixing operating below the convective envel ope durin g the RGB 
phase (the so called extra-mixin}^, see e. g. iNoUett et al.. ,20031 : 
ICharbonnel & Zahnll2007t iDenissenkov & Pinsonneault I2OO8I) . 
Our reference model (ST, dash-dotted magenta curve) has an 
initial mass M = 1 .5 Mq with Z = 0.006, corresponding to 
[Fe/H] = -0.36. We assume an initial solar-scaled composition, 
which implies C/O = 0.5 and '^C/'^C = 90. After the FDU, 
the surface values attain C/O = 0.36 and '^c/'^c = 24. This 
is due to the fact that the convective envelope penetrates into re- 
gions where partial hydrogen-burning occurred before. No extra- 
mixing has been included in the ST case. The values so obtained 
clearly disagree with the M stars observations, which show an 
average C/O = 0.18 and an average '^C/'^C - 11. Thus, we 
explored the possibility of an occurrence of extra-mixing on the 
RGB. 

This hypothesis is supported by the bulk of observations of 
RGB stars in the galactic field, as w ell as in open and globular 
clusters (see e. g. lGratton et al.ll2000l) . These observations show 
that this additional mixing occurs in low-mass stars (M < 2 Mq) 
during the first ascen t along the Red Giant Branch. Moreover, 
lEggleton et al.l ( l2008l) identified a mixing mechanism driven by 
a molecular weight inversion (Sfi-mixing) in three-dimensional 
stellar models that must operate in all low-mass stars while they 



' The errors for (J - K) are higher, ranging up to and being even 
greater than 0. 1 . This could also affect the differences between the ef- 
fective temperature scales discussed in Sect.|4] 
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Fig. 7. Comparison between observational data and theoretical models. '^C/'^C isotopic ratios vs. C/O ratios are reported. In the 
left panel we report the range of O-rich stars only, while in the right panel the whole sample has been considered. The star symbols 
with error bars indicate our observations. The red dotted, the blue dashed, and the green long-dashed curves in the left panel are 
partly covered by other tracks in the plot. The symbols in the left panel mark the end point of the tracks after the RGB phase. In the 
right panel, the symbols indicate the values attained during the interpulse phases, which are on the order of 10^ times longer than 
the TDU episodes. See text for details. 



are on the RGB. The operation of this RGB extra-mixing is also 
required to explain the relatively low '^C/'^C ratios in the M 
stars of NGC 1846 (paper I). 

As in paper I, we include this additional mixing right after 
the RGB luminosity bump. The extension of the zone in which 
this additional mixing takes place is fixed by prescribing the 
maximum temperature the material is exposed to (T'^gb-'- 
circulation rate is tuned by setting the mixing velocity to a value 
that is a small fraction (cf. paper I) of the typical convective ve- 
locities in the envelope of an RGB star As discussed in paper I, 
the carbon isotopic ratio largely depends on T^^^. In Fig. [T] we 
report two models, characterised by T^^'^ = 2.0 x 10^ K (red 
dotted curve) and = 2.5 x 10^ K (blue short-dashed curve). 
In the first case, the final '^C/'^C ratio decreases, reaching a 
value in good agreement with those observed in the M stars of 
NGC 1978. Notwithstanding, the C/O ratio remains unaltered 
and higher than the observed one. An increase of T'^qb up to 
2.5 X 10^ K leads to a lower C/O ratio of about 0.21, which is 
in better agreement with the observations. However, the corre- 
sponding '^C/'^C ratio is 5, which is definitely lower than the 
average value (11). Note that only LE5 shows such a low value 
(i^C/i^C = 6). 

Then, we explored the possible effect of an oxygen enhance- 
ment (see the discussion in Sect. 12. It . The black solid line 
in Fig. |7] represents a model similar to the ST case, but with 
[O/Fe] = 0.2. As for the ST case, the effects of the FDU are 
clearly recognisable. An additional model, as obtained by in- 
cluding an RGB extra-mixing (r^cg = 2.1 x lO^K) is also re- 
ported (green line). The final surface composition of this model 
(C/O = 0.23 and '^C/'^C = 13) is close to the average val- 
ues of the observed sample (C/O = 0.18 and '^c/'^c = 11). 
Thus, a first conclusion is that moderate RGB extra-mixing 



(2.1 X 10^ < r™^ < 2.3 X 10^ K) coupled to moderate oxygen 
enhancement appears to nicely reproduce the observed average 
composition of the M stars in NGC 1978. The same conclusion 
has been reached in the case of NGC 1846 (paper I). 

What are the implications of this scenario when applied to 
the whole observational sample (O-rich and C-rich stars)? We 
carry on our analysis assuming that the carbon stars are intrinsic, 
i. e. that their surface composition, in particular the high C/O 
and carbon isotopic ratios, is the result of nucleosynthesis and 
mixing processes occurring during the thermally pulsing AGB 
evolution. 

In the right panel of Fig. [T] we report the same models as 
shown in the left panel. None of the theoretical tracks can si- 
multaneously reproduce the abundance ratios of the M and C 
stars. The high values of the C-star carbon isotopic ratio can 
only be reproduced by the oxygen-enhanced model with no RGB 
extra-mixing. All the other models predict too low '^C/'-'C ra- 
tios. The situation appears even more peculiar when compared 
with the case of NGC 1846, for which we found C stars with 
higher values of C/O, but with ^^C/^^C between 60 and 70, 
so that an AGB extra-mixing (in addition to the RGB extra- 
mixing) was invoked to reproduce the observations. The vast 
difference between the two clusters is illustrated in Fig. [8] In 
the left panel, data for NGC 1978 are compared with three mod- 
els computed under different assumptions for the RGB and the 
AGB extra-mixing, namely no extra-mixing (black-solid line), 
moderate RGB extra-mixing only (green-long-dashed line), and 
moderate AGB extra-mixing only (blue-dashed line). All the 
three models have M - 1.5 Mq and [O/Fe] = 0.2. In the right 
panel, data for NGC 1846 are compared with similar models. 
In this case, the mass is 1.8 and the red-dashed line refers 
to a model with both RGB and AGB extra-mixing (see paper I 
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for more details on models with AGB extra-mixing). In sum- 
mary, the high carbon isotopic ratio observed for the C stars of 
NGC 1978 ('^C/'^C > 150) rules out both RGB and AGB extra- 
mixing, which are instead required to reproduce the evolutionary 
sequence of M, S, and C stars in NGC 1846. 

5.2. Alternative scenarios 

In the previous section, we did not identify a theoretical sce- 
nario that satisfactorily reproduces the whole set of stars be- 
longing to NGC 1978 (including both the O-rich and the C- 
rich sub-sample). Moreover, the intermediate-age LMC cluster 
NGC 1846, which we investigated in paper I, presents a totally 
different abundance pattern, which cannot be ascribed to the 
same theoretical scenario. The lack of S stars is an additional 
peculiarity of NGC 1978. So far, our arguments have been based 
on the assumption that the stars in our sample possess the same 
initial composition, mass, and age so that they form an evolution- 
ary sequence, in which the variations of the surface composition 
are due to intrinsic processes only. In this section, we speculate 
about possible alternative scenarios, namely: 

1 . the existence of two stellar populations having different ages; 

2. the existence of rejuvenated stars; 

3. the activation of additional mixing at the base of the convec- 
tive zones generated by a thermal pulse. 

Th e existence of two populations proposed bv lAlcamo et al.l 
(1199 9') is an attractive possibility, because it might explain the 
dichotomy between the M and the C stars. The M stars would 
then belong to the older population. Their mass is sufficiently 
low for the occuiTence of the RGB extra-mixing, but too low for 
the occurrence of the TDU during the AGB. On the contrary, the 
C stars are younger (less than 1 Gyr) and more massive (M > 
2 Mq). For these stars, the persistence of the /i-barrier at the base 
of the hydrogen-rich envelope prevented the RGB extra-mixing, 
whilst, during the AGB, TDU accounts for the observed large 
C enhancement. However, this scenario seems unlikely. First of 
all, a number of other studies did not find significant metallicity, 
mass or age differences. Second, a cluster merger creating two 
populations is not co nsidered to be a valid option (see Sect. |2]i. 
lAlcafno et af] (Il999 l) argue that the north-western half could be 
more metal-rich and maybe slightly younger. The distribution of 
our targets in Fig. [T]does not support this idea. The M stars are 
distributed equally around the cluster centre. The carbon stars in 
our sample lie either in th e centre or in the south-eastern half of 
the cluster. According to iLlovd Evans! ([19801) . there is another 
carbon star (LEll) located in the north-western half. Although 
the C stars appear to line up along the cluster's major axis, there 
is no indication that they concentrate in a certain region. A wide 
spread in metallicity would have shown up in the fitting process, 
yet we were able to obtain reasonable results for both oxygen- 
rich stars and the hottest carbon-stars with a single value for the 
metallicity. 

The scenario of a single age stellar population does not ex- 
clude the existence of rejuvenated stars. Merging, coalescence, 
or mass accretion processes are often invoked to explain the Blue 
Stragglers observed in many GCs. Actually, a non-negligible 
number of stars lie above the turn-off and the sub-giant br anch in 
the colour-magnitude diagram of NGC 1978 (.MucciarelU et alj 
I2007ah . although we cannot exclude that those are field stars 
or the result of unresolved pairs. If the product of one of the 
above mentioned processes is a star with a mass larger than 
2Mo, or if the rejuvenating process leaves the envelope com- 
position enriched with heavy elements, the efficiency of the 
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Fig. 9. '^C/'-'C isotopic ratio vs. C/O ratio for a series of mod- 
els including an overshoot at the base of convective zone gener- 
ated by the TPs. The extension (in fraction of the pressure scale 
height) of the overshoot zone is indicated in the legend. 

RGB extra-mixing is significantly reduced (ICharbonnel & ZahnI 
12007 1). Such a scenario might explain the high carbon isotopic 
ratio of the C stars, but the lack of S stars remains a mysteryQ 
Moreover, even if we do not have a detailed statistics of the num- 
ber of binary systems in NGC 197fl it seems unlikely that two 
merger events occur nearly simultaneously to produce the stars 
B and LE6, for instance. 

The last alternative scenario we have investigated regards the 
possible modification of the efficiency of the mixing processes 
taking place during the AGB phase, which is responsible for the 
increase of '^C and in the envelope. A deeper TDU would 
increase both the carbon isotopic ratio and the C/O ratio, thus 
leaving the slope of the curves in the right panel of Fig. |7]unal- 
tered. In this scenario, the possibility of producing S stars as a 
consequence of the TDU is reduced, since the transition to the 
C-star stage is more rapid. We recall that the TDU depth ba- 
sically depends on the strength of the thermal pulses: stronger 
TPs push outward the external layers more easily, therefore in- 
ducing deeper dredge-up episodes. In turn, the strength of a ther- 
mal pulse depends on many physical and chemical properties, 
such as the c ore and the envelope m ass, or the CNO content in 
the envelope dStraniero et al.ll2003b . On the other hand, an in- 
crease of the primary oxygen that is dredged up might limit the 
increase of the C/O at the surface and steepen the '^C/'^C-C/O 
relation. The fraction of oxygen produced by the '^C + a reac- 
tion increases toward the centre. Therefore, an additional mixing 
process, as due to a ballistic overshoot or to non-standard mech- 
anisms (such as rotation or magnetic kick), that moves primary 
O (and C) from the deep interior to the He-rich inter-shell, may 



Note that our sample contains almost all AGB stars in this cluster 
(see Fig.|6j. 

' Estimated binary fr actions in other LMC and SMC clusters range 
from 10 (NGC 1783. iMucciarelli et al.l l2007bh to < 30 per cent 
jChiosi & Vallenarill2007l) ! 
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Fig. 8. Comparison between observational data and theoretical models. '^C/'^C isotopic ratios vs. C/O ratios are reported. In the 
left panel, we report observational data of the cluster NGC 1978, accompanied by selected theoretical tracks. As a comparison, in 
the right panel, we report the fitting curve for another LMC cluster (NGC 1846). The symbols along the tracks indicate the values 
attained during the interpulse phases (see also Fig.|71l. The stai" symbols with error bars indicate our observations. See text for details. 



indirectly explain the high carbon isotopic ratio of the carbon 
stars in NGC 1978. To shed light on this possibility, we com- 
puted some additional models, in which we applied an artificial 
overshoot at the base of the convective zone generated by the 
TPs. The initial mass adopted in the models is 1.5 Mq, and the 
initial composition is oxygen-enhanced (-1-0.2 dex). A moderate 
RGB extra-mixing (T'^q^ - 2.1 x 10^ K, corresponding to the 
case of the green long-dashed line in Fig.|7]i has been included. 
We parametrise the additional mixing by limiting the extension 
of the overshoot layer to a fraction of the pressure scale height 
of the most internal convective mesh point. In Fig. |9] we report 
the resulting '^C/'^'C isotopic ratios as a function of C/O. 

Even if these models apparently provide a simultaneous re- 
production of the composition of both M and C stars, we have to 
stress the theoretical and observational evidence against such a 
process. At the He-flash peak, the inner border of the convective 
zone already reaches the layer where the temperature attains its 
maximum value. Below this point, due to the neutrino energy 
loss, the temperature (and the entropy) decrease s, so that the 
decele ration due to the buoyancy is quite strong. iHerwig et alJ 
(12006 1) showed that even if the pressure and the entropy barrier 
at the base of the convective shell generated by a TP is stiff, con- 
vective plumes can penetrate into the underlying radiative layers 
due to g-mode oscillations. These structures reach down to about 
2x10^ cm below the formal convective border (see their Fig. 24), 
and develop mean vertical velocities that are about 1000 times 
lower than the average convective velocities. The extension of 
this penetration has the same order of magnitude as the over- 
shoot zone in our models, namely 1.0 x 10^ cm, 1.2 x 10^ cm, 
and 1.5 x 10^ cm, for 0.25 Hp, 0.35 Hp, and 0.5 Hp, respec- 
tively (these are mean values). However, if in our ID models 
the whole mass of the spherical overshoot region is efficiently 
mixed, the amount of mass contained in the penetrating plumes 
of the 3D model is certainly lower. In addition, the result of the 



3D model may be a ffected by the assum ed boundary conditions. 
As was stressed bv iHerwig et al.l ( 12006 ). their models ignore fi- 
gradients, which increase the stability and reduce mixing across 
the formal convective boundaries. Moreover, they found that the 
amplitude of the gravity waves depends both on the resolution 
adopted in the simulation and on the choice of the heating rate. 

On the other hand, the bulk of the C stars in the Milky Way 
as well as those in NGC 1846 have '^C/'^C rati os in the range 
from 40 to 70, and only in a few cases up to 100 dLambert et al.l 
il986h . This occurrence leads us to the conclusion that such an 
overshoot is supposedly an uncommon process. Note that the 
occurrence of an overshoot below the convective zone generated 
by a TP causes considerable changes of the physical conditions 
in the He-rich inter-shell. In particular, higher temperatures de- 
velop at the pulse peak and, therefore, the ^^Ne(a, n)^^Mg reac- 
tion rate becomes an important neutron source, even in low-mass 
AGB stars. In such a case, the resulting i-process nucleosynthe- 
sis would be characterised by an overproduction of neutron-rich 
isotopes. For instance, the expected isotopic compositions of 
strontium, zirconium, molybdenum, and barium would substan- 
tially differ from those measured in pre-solar SiC grains, which 
were produced in the cool circumstellar en velope of a past gen- 
eration of C stars (see iLugaro et alj[2003[) . Additionally, calcu- 
lated element abundances ratios such as Rb/S r would not matc h 
observations obtained from carbon stars ('see lAbia et al.ll200ll) . 
For these reasons, if the overshoot from the base of the convec- 
tive zone generated by TPs is non-negligible, it should be a rare 
event. 

In summary, a clear and coherent picture of the AGB stars 
in the cluster NGC 1978 within the more general context of low- 
mass AGB stellar evolution and nucleosynthesis cannot be eas- 
ily achieved. Once again, we stress the uniqueness of the C star 
sample of NGC 1978 with respect to other LMC clusters. 
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6. Conclusions 

In this paper, we presented a sample of AGB stars belonging to 
the LMC cluster NGC 1978. Where it was possible, we derived 
the C/O and the '^C/'^'C ratios. The spectroscopic data reveal 
the presence of two sub-samples, one where C/O < 1 in the stel- 
lar atmosphere (oxygen-rich or M-type stars), and the other with 
C/O > 1 (carbon-rich or C-type stars). The oxygen-rich stars 
present low values for the C/O ratio (mean value 0.18) and the 
'^C/'^C ratio (mean value 1 1). The spread in the C/O ratio is re- 
stricted, all values are consistent with the mean value within the 
error bars. We observe a wider spread in the '^C/'-'C ratios. We 
conclude from the spectroscopic data in comparison to our evo- 
lutionary models that the TDU mechanism is not (yet) working 
in these stars. The carbon-rich stars of our sample instead show 
C/O values and '^C/'^C ratios consistent with the occurrence of 
TDU. A lack of S stars, with C/O ratios between those of the 
O-rich stars and the C-rich stars, has been highlighted. 

We did not find a theoretical scheme that is able to satisfacto- 
rily reproduce the chemical abundance pattern in NGC 1978. By 
claiming the existence of a non-standard mixing mechanism, we 
postulated some possible solutions. However, we are aware that 
too many ad hoc assumptions make our analysis objectionable, 
taking into account that some of them are strongly limited by 
observational constraints. In particular, a more consistent picture 
of all the additional mixing mechanisms, which are active during 
the RGB and the AGB phase, is required. We also discussed the 
scenario of multiple stellar populations and rejuvenated stars, 
and highlighted the differences to the abundance pattern in the 
cluster NGC 1846. 

Narrower constraints for stellar evolutionary models could 
be derived from atmospheric models including a set of accu- 
rate molecular line data, especially in the case of carbon stars. 
Additionally, the more evolved a star is (thus having a higher lu- 
minosity and C/O ratio), the more dynamic effects influence the 
spectral appearance. From the observational point of view, im- 
provements could be achieved by means of the choice of the ob- 
served wavelength regions, some of which are less crowded with 
molecular lines. Moreover a comparison between low-resolution 
data and high-resolution data should be explored further (cf. 
IWahlin et al. 2006). If low -resolution data deliver robust C/O ra- 
tios (and stellar parameters), the isotopic ratios could be derived 
in the high-resolution spectra more accurately. 

In conclusion, we affirm that there is much room for im- 
provement. As a first step, we are currently acquiring new data 
to increase our statistics by studying more AGB stars in other 
LMC globular clusters. In this way, we could shed light on the 
puzzling abundance pattern that we found in NGC 1978. 
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Appendix A: Changes to the CO line list 

To match the positions of the CO lines in the recorded K-band 
spectra, we had to modify s ome transition frequencies as given in 
the list of lGoorvitchI (|1994|FI The wavenum bers were taken from 
the Infrared Atlas of the Arcturus Specti-um jHinkle et al.lll995l) . 
The changes are documented in Table lA.l"! The given quantities 
are the old and new transition frequency (cr), the lower state term 
energy E\, the logarithm of the value (\og gf), the lower and 
upper state vibrational quantum numbers (vi and Vu), the lower 
state angular momentum number { J\), and the isotope identifica- 
tion. 

At a later point in the analysis, we compar ed the correc- 
tions of the line positions to the CO data given bv lChandra et alj 
([1996). The wavenumbers from the Arcturus atlas are in good 
agreement with their transition frequencies. The differences are 
less than 0.04 cm"' for the 1-3 vibration transitions and less 
tha n 0.01 cm"' f or the other lines. However, we did not use 
the IChandra et al.l ( 1996) data, because the line strengths of the 
'^CO 3-0 band head in the H band from this list are severely 
underestimated and therefore unusable for the analysis. 
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Table A.l. Modified line positions of CO lines in the K band 



r -In 
cr„M [cm '] 


r -In 
O-new [cm '] 


£1 [cm"'] 




Vl 




-/i 


Isotope 


4219.5583 


4225.3630 


19009.6550 


4.865E-05 


1 


3 


95 


r ! ^ 1 f, 


4226.9718 


4227.7190 


14547.6172 


6.569E-05 


2 


4 


74 


1 T /-^ 1 ft/^ 

'-C'O 


4223.6261 


4228.9450 


18668.6535 


4.782E-05 


1 


3 


94 


12/-- lfi/"v 


4228.9365 


4229.6010 


14278.1804 


6.439E-05 


2 


4 


73 


12/-- lfi/-v 


4230.8124 


4231.4030 


14012.1198 


6.310E-05 


2 


4 


72 


12/-- !fi/-v 

''^C 


4227.5739 


4232.4430 


18330.8030 


4.701E-05 


1 


3 


93 


'^C 


4232.1241 


4233.3140 


11962.5127 


1.167E-05 
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81 




4234.3009 


4234.7650 


13490.1716 


6.057E-05 


2 


4 


70 


12^160 


4234.3757 


4235.4470 


11676.5472 


1.144E-05 
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80 


13^160 


4235.9151 


4236.3250 


13234.3055 


5.932E-05 


2 


4 


69 


12^160 


4237.4434 


4237.8040 


12981.8588 


5.809E-05 
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68 


12^160 


4238.8865 


4239.2040 


12732.8418 


5.687E-05 


2 
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67 


12(-.16o 


4240.2450 


4240.5230 


12487.2648 


5.566E-05 
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66 


12(-i16q 
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4241.7630 


12245.1376 


5.447E-05 
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65 
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